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Abstract. We show, in general, how to transform the nonautonomous nonlinear Schrodinger equa- 
tion with quadratic Hamiltonians into the standard autonomous form that is completely integrablc 
by the familiar inverse scattering method in nonlinear science. Derivation of the corresponding 
equivalent nonisospectral Lax pair is also outlined. A few simple integrable systems are discussed. 



1. Introduction 

Recently several nonautonomous (with time-dependent coefficients) and inhomogeneous (with 
space-dependent coefficients) nonlinear Schrodinger equations have been discussed as (possible) 
new integrable systems [7J, [8J, [H], [IB], [20], 02], 08], [52], [53], [59], [66], [67], [76], [8T], [82], 
[83], [H], [85], [87], [93], [96], [97J, [98], [M] (see also [2], [3], [9], [18], [19], [21], 05], [55], [92] and 
references therein for earlier works). They arise in the theory of Bose-Einstein condensation [35J, 
[75] . fiber optics [6], 07J, superconductivity and plasma physics [18], [19], [70], [TTj . 

As pointed out in recent papers [7j, [IB] . 08] . 09] and [55] (see also [2], [3], [21], [31], 05] . 
[68] , |74j ) . these systems can be reduced by a set of transformations to the standard autonomous 
nonlinear Schrodinger equation, which explains their integrability properties because this equation 
is a well-known complete integrable system with Lax pair [38], |101j . |102j . [103] , conservation 
laws and A^-soliton solutions, solvable through the inverse scattering method [2], [3], [5J, [33], [71] . 
[78] . [101] . |102j . |103j . Integration techniques of the nonlinear Schrodinger equation include also 
Painleve analysis [13], [23], [23J, [2S], [21], [33], [31], 09], [53], jSS], [89], [93], Hirota method 
[31] . [31], Backlund transformation [12], [T7J, [M] and Hamiltonian approach PQ, [3J, 03], 03] - 
[63], [71] among others [36], [62], [72J, [TTJ. 

A goal of this paper is to construct these transformations explicitly (in quadratures) for the 
most general variable quadratic Hamiltonian. A simple relation with Green's function of the linear 
problem, which seems has been missing in the available literature, is emphasized. Basics of the 
classical soliton theory including one and two soliton solutions, inverse scattering technique and the 
corresponding equivalent Lax pair are also briefly summarized in order to make our presentation as 
self-contained as possible. It may facilitate applications of our result to specific nonlinear integrable 
systems. 
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2. Transformation into Autonomous Form 



We consider the nonautonomous nonlinear Schrodinger equation: 

t^ = H^ + h\^\ 2 ^ (2.1) 

on K, where the variable Hamiltonian H = Q (p,x,t) is an arbitrary quadratic form of operators 
p = —id/dx and x, namely, 

iif> t = -a (t) i) xx + b (t) x 2 i, - ic (t) xiJ x -id{t)iP-f (t) x^j + ig (t) + h (t) |V>| 2 V, (2-2) 

(a, b, c, d, f and g are suitable real- valued functions of time only) under the following integrability 
condition^] [87]: 

h = h a (t) (3 2 (t) ii it) = h (3 2 (0) p 2 (0) a[t)> v it) (2.3) 

p [t) 

(ho is a real constant, functions (3, A and p will be defined below). 

We present the following result. 
Lemma 1. The substitution 

1> (x, t) = _A^ e K^ + s(t )x+K (t)) x ^ r) j z = p(t) x + e (t), r = 7 (t) (2.4) 

transforms the nonautonomous and inhomogeneous nonlinear Schrodinger equation 112. || ) into the 
standard autonomous form with respect to new variables £ = (3 (t) x + e (t) and r = 7 (t) : 

«Xr + ^0 \x\ 2 X = Xtf (2-5) 



provided that 



and 



where 



^ + 6 + 2ca + Aaa 2 = 0, (2.6) 

^ + (c + 4aa)/3 = 0, (2.7) 
dt 

^ + a(3 2 = (2.8) 
at 

— + (c + 4aa) 0" = / + 2ac/, (2.9) 

(JjL 

de 

- = ( g -2a6)P, (2.10) 

^=^-a5 2 , (2.11) 

4a (t) p (t) 2a (t) v ; 



1 If the nonlinear term has the form h \ip\ p tp, popular in the mathematical literature, the condition becomes 
h = h a ap 2 pPl 2 . 
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The autonomous equation (12.51) is completely integrable by advanced methods of the soliton 
theory [5J, [53], [7_T], [79], |1U1] . |102j . |1U3] (see also [36J and references cited in the introduction). 
Equations (I2.2p -( l2~3l) seem to represent the maximum nonautonomous and inhomogeneous one- 
dimensional integrable system of this kind. (Important special cases of the transformation (12 .4p are 
discussed in Refs. [2], [15], [19], [2TJ, 05], @H], [55], [73], [76] and [87].) 

Our transformation (12.41) involves the real-valued functions a, 0, 7, 5, e and k (of time t only) 
defined as solutions of the system of ordinary differential equations (I2.6p - (I2.11I) . This nonlinear 
ODE system has already appeared in Ref. [28] from a different perspective (we shall refer to this 
system as a Riccati-type system). The substitution (12.121) reduces the Riccati equation (12.61) to the 
second order linear equation [25] : 

H 1 ' - r (t) 11' + Aa (t) [i = (2.13) 

with 



d 



r(t) = 2c + Ad, a (t) = oh - cd + d 2 + 

" d 



d' 



( 2 - 14 ) 

a 2 
(Relations with the corresponding Ehrenfest theorem for the linear Hamiltonian are discussed in 
Ref. [30] . It provides a clear physical interpretation of our results in case of the Gross-Pitaevskii 
model of Bose condensation.) Equation (12.81) implies the monotonicity of the new time variable 
r = 7 (t). 



Proof. Differentiate i/j = /i~ 1/2 (t) e iS{x < t] x (f , r) with 5 = a (t) x 2 + 5 (i) x + k (t) and ^ 
£(*), r = 7 (f) : 



1 



- («x 2 + 5'x + k') x + i ( {fix + e') Xi + l'Xr ~ 



-iS 



[2 (2as + 5) x + PX(] 



and 



-i.S 



4= [( 2 ^« - ( 2 « s + 5 ) 2 ) X + 22 (2ax + 5) ^ + f?Xt J 
v A* 



/? (*) x + 

(2.15) 
(2.16) 

(2.17) 



Substitution into (12. 2p . with the help of the integrability condition (12. 3 p and the system (I2.6j) ( T2TTTj) . 
results in (12. 5p . Computational details are left to the reader. □ 

This observation provides a new interpretation of the Riccati-type system (I2.6j) f T2TTT]) . which 
has been originally derived in Ref. [28] during integration of the corresponding linear equation via 
Green's function. 



3. Integration of the Riccati-type System 



In order to construct the transformation (12.41) explicitly, one has to solve the nonlinear ODE 
system (I2.6p - (l2.1ip . As already known, the initial value problem of the Riccati-type system, which 
corresponds to the linear Schrodinger equation with a variable quadratic Hamiltonian (generalized 
harmonic oscillators [H], [40], [46], [95], [100] ). can be explicitly solved in terms of solutions of 
characteristic equation (12 . 131) as follows [28], [HD], [HE], [88] : 

fi (t) =2fi(0)fi (t) (a (0) + 7oW), (3.1) 
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and 



Here, 



and 



a (t) = a (t) 
(3(t) = 



4(«(0)+ 7o (*))' 
P{0)/3 (t) _ /3(0)/x(0) 
H (t) 



t(*) = t(o) 



2 (a (0)+ 7o (*)) 
/3 2 (0) 



A(t), 



K{t) 



So (t) 

6(0) 



A (a (0) + 7o (t)) 

go (t) (o(0)+ go (t)) 

2(a(0)+7o(*)) ' 
/3(0) (g(0)+ £o (t)) 

2(a(0)+ 7o (t)) ' 

(<f(0)+M*)) 2 



« (0) + « (t) 



4(a(0)+ 7o (*))' 



a (t) 
7o (*) 



1 ^(t) d(t) 



4a (t) fjtQ (t) 2a (t) 
X(t) 



fi (t) 
1 



A (t) = exp 
Mi(0 + d(0) 



c s 



2d(s)) ds 



2M0)// (t) 2a (0) 



5 (*) 



A(f) 



Mo (*) 



L 



f(s)-TTZ9(s))n u {s) - 



a is 



2a (s) 



Mo ( s ) 



ds 

W) 



2a (t) A (t) 



5o (*) + 



+2 



Mo (*) 

4 a(s)A(s) 



Mo ( s ) 



a (s) a (s) A (s) 

o (Mo( s )) 2 
N d(s) 



a s 



O (t) Mo (*) r2 



Mo (*) 



*S (0 - 4 



a (s) a (s) 
o (Mo ( s )) 2 



(/io (s) <J (s)) ds 
ds, 

(/i (s)o" (s)) 2 ds 



-2 



as 



(^ (s)5 (*)) /( 



d(s) 



as 



ds 



(3.2 
(3.3 
(3.4 

(3.5 
(3.6 
(3.7 

(3.8 
(3.9 
(3.10 

(3.11 
(3.12 



(3.13) 



'o Mo( s ) 

(So (0) = — e (0) = g (0) / (2a (0)) and k (0) = 0) provided that fx Q and n x are the standard solutions 
of equation (I2.13P corresponding to the following initial conditions /x (0) = 0, /i (0) = 2a (0) ^ 
and fjL t (0) 7^ 0, fi[ (0) = (proofs are outlined in Refs. [28], [32] and [HS]). (Formulas (l3T8l)- fl37[3|) 
correspond to Green's function of generalized harmonic oscillators; see, for example, [25], [SO], |41j . 
[56] , [5"U] . [SB] , [55] and references therein for more details.) 

One may refer to the solution fl3.8p - fl3.13p as the fundamental solution of Riccati-type system 
fl2.6p - fl2.lip . Thus the transformation property fl3.lD - fl3.7p allows one to find solution of the initial 
value problem in terms of the fundamental solution (for the nonlinear ODE system). 
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4. Integration of the Nonautonomous Linear System 

The transformation (12. 4p reduces the linear Schrodinger equation of generalized harmonic oscil- 
lators, namely, equation (12.11) with h — 0, to the Schrodinger equation for a free particle i\ T = 
with a familiar Green function given by 

\2 



1 



. exp 



'4(r-r ) 



(4.1) 



y-Airi (r - r ) 

where £ = (3 (t) x + e (t) , r\ = (3 (0) x + e (0) and r = 7 (£) , r = 7 (0) . One can verify directly that 
Green's functions of generalized harmonic oscillators [28] . 

G (x, y, t) = z exp [i (a (t) x 2 + /3 (t) + 7o (t) + 5 (t) x + e (t)y + k q (t))] , (4.2) 
V 27r ^o W 

are derived from the simplest free particle propagator (14. ip with the help of our transformations 
(I3.ip - (l3.7p . In is worth noting, though, that the transformation (12. 4 p requires a knowledge of the 
functions /x, a, (3, 7, 5, e and k, which allows to determine the Green's function for the generalized 
harmonic oscillators directly from (12. 2p . Thus finding of this transformation is equivalent to inte- 
gration of the original linear equation from very beginning. (Lemma 1 extends this observation to 
the nonlinear Schrodinger equation ( j2.2h ) . 

Then the superposition principle allows to solve the corresponding Cauchy initial value problem: 



for suitable initial data ip ( x , 0) = (f (x) (see Refs. [28], [EE] 
As shown in [86J, the following asymptotics hold 

1 c(0) a'(0) 



G(x,y,t)ip(y,0) dy 

] and [86J for details). 



(4.3) 



a (t) 

7o 00 
So (t) 



4a (0) t 4a (0) 8a 2 (0) 



+ 0(t) 



(4.4) 



a' (0) 

2a(0)t ' 4a 2 (0) 



+ 



0(f), 



+ 



c(0) a'(0) 



4a (0) t 4a (0) 
9(0) 



2a (0) 
«o (*) = O (t) 



+ 0(t) 



8a 2 (0) 



+ 0(t), 

9(0) 
2a (0) 



+ o(t) 



as t — > for sufficiently smooth coefficients. Then 



G(x,y,t) 



1 



'27170 (0) t 



x exp 



exp 



. ip -yf 
% 4a (0) t 



(4.5) 



a' (0) 



(x-y) + 



c(0) 



(x 2 - y 2 ) 



9(0) 



(x-y) 



^a 2 (0) v y/ 4a (0) v y y 2a (0) 
as t — y 0, which corrects a typo in Ref. [28]. (Here, / ~ g as £ — )■ 0, if lim^o (f 1 9) = !•) 
Another form of solution of the linear problem can be found by an eigenfunction expansion 
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Numerous examples of (super) integrable (driven) generalized harmonic oscillators with a de- 
tailed bibliography can be found, for instance, in recent publications [23J, [2H], [3D], [HE] and [32J. 
In addition, our Lemma 1 shows that these results provide explicit transformations (12.41) of the 
corresponding nonlinear systems into the standard completely integrable forms. 



5. One Soliton Solution 

In the next few sections, we summarize basics of the inverse scattering technique for the au- 
tonomous nonlinear Schrodinger equation ( 12. 5 ft in order to facilitate use of the transformation ( 12.41 ) 
for specific nonautonomous and inhomogeneous nonlinear Schrodinger equations ( 12. 2p from various 
applications. More details, when needed, can be found in classical works [2J, [3J, [33], [71], [79], 
[101] ■ [102] . [103] (see also references cited in the introduction). 

As well-known, equation (12. 5ft has a travelling wave solution of the form 

X (£, r) = e <W^o)+0) F (e + 2ry) (51) 

provided 

(dF\ 2 1 
— — I = Co + goF 2 H — h F A (Co is a constant of integration) . (5.2) 
dz J 2 

Examples include bright and dark solitons, and Jacobi elliptic transcendental solutions for solitary 
wave profiles [2J, [33], [TTJ, [79], [STJ. Setting C = y = 0, gives the stationary breather, which is 
located about £ = and oscillates at a frequency equal to g [75], [79] . 

By (12.41) . the nonautonomous Schrodinger equation ( 12. 2ft under the integrability condition (12. 3ft 
has the following solution: 

■0 (x,t) = — — exp (i (ax 2 + + 7 (y 2 — g Q ) + 5x + ey + «)) (5.3) 
v ^ 

xF(/3x + 27?/ + £), 

where the elliptic function F satisfies equation ( 15.2ft and 0, g>Q and /to are real parameters (see 
also Ref. [87] for a direct derivation of this solution). 

6. Integrability of Nonautonomous Nonlinear Schrodinger Equation 

The substitution (X, T) = v^oX (\%X, — 2T\ transforms equation (12.51) into the standard 
forms 

m> T + m xx ±2\m\ 2 m = o (6.1) 

(focusing and defocusing), which can be obtained as the flatness condition: 

Ut-Vx + UV -VU = (6.2) 
for the Lax-(Zakharov-Shabat) pair: 

U = -iA<7 3 + * <r+ T (6.3) 
-iA * 

and 

V = i(-2\ 2 ±\ty\ 2 )a 3 + (2\ty + ity x )a + ±(-2\ty* +ity* x )a_ (6.4) 
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' 






f -i \ 




"' = 1 




o • 


0"2 = | 


v i ' 


^3 = f 



i(-2A 2 ±|^| 2 ) 2A^ + ^ X 
T2\y*±M* x i(2A 2 ^|#| 2 ) 

(we use the asterisk for complex conjugation). Here, A is a constant, a± = {a\ ± ia 2 ) /2 and cti, 
<j 2 , cr 3 are the Pauli matrices: 

(6 - 5) 

Since the Lax pair guarantees complete integrability and can alone derive all its associated prop- 
erties, our transformation fl2.4jl trivially explains the integrability features of the nonautonomous 
nonlinear Schrodinger equation (12.11) . including X-soliton solutions, infinite conservative properties, 
etc., (see Refs. [53], [7T] and [7H] for more details). 

7. Inverse Scattering Method 

Solution of the Cauchy initial value problem through the inverse scattering method is discussed 
in [3], [5], P], HQ, [79], HDI], HD2] and [TU3]. In the focusing case, 

^ T + ^ xx + 2|^| 2 * = 0, (7.1) 

the Zakharov-Shabat system contains four equations for an auxiliary two-component wave function 
$ = (<p,vf : 

#x = # T = V$, (7.2) 

namely, 

cp x = -iXip + ^v, (7.3) 
vx = —^*<p + iXv 

and 

y? T = i(-2A 2 + |^| 2 )y? + (2A^ + ^ x )^, (7.4) 
v T = (-2A^* + M* x )<p + i (2A 2 - |^| 2 )^. 
Assuming that * (X, T) — >• (and ^ x (X, T) -)> 0) and as X ->> ±oo implies 

yj T -2i\ 2 (p, i; T -»> 2zA 2 i; (X -> ±oo) (7.5) 
so the scattering data for the problem 

i(^=A(rV (7-6) 



VI \ V 



L = i<T S -^-i*a + -iVa- = i(% _*) (7.7) 

evolve with time as 

6 (A, T) = 6 (A, T ) e^Cr-ib), r „ ( T ) = r„ (T ) e «*<r-H (7.8) 

Then Cauchy initial value problem for the nonlinear Schrodinger equation (16. ip can be solved as 
follows [5J], [S], Uni], HD2], [103] : 

*(X,T) = -2X(X,X,T), (7.9) 
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where K (X, Y, T) satisfies the linear integral equation 

K(X,Y,T) = B*(X + Y,T) (7.10) 

POO POO 

- / K (X, Z, T) B (Z + W, T) B* (Y + W, T) dZdW 
J x Jx 

and B (X, T) can be obtained in terms of the scattering data: 

N -i /*00 

B (X, T) = -zJ2 r n (T ) e <^x+^(T-T )) + J_ b (A> Tq) ^xx+rtp-n,)) dX (7 n) 

n=l 271 

(see Refs. 0, [3], [5], [M], [H], [ZH], [79], [91], [M], [102], [103] for more details). 
As a result, the following combination of gauge, scaling and coordinate transformations [55] : 

V> *) = 7^y e<a(i)x2+5(<)a;+KW) 117 (71 ^ W x + £ ' ~h {t) ) ' (7 - 12) 

written here explicitly in terms of solution of the Riccati-type system (I2.6p -( 12TT|) from section 3, 
allows to solve Cauchy initial value problem for the nonautonomous nonlinear Schrodinger equation 
(12. 2p with the help of the standard inverse scattering technique. The following choice of parameters: 
a (0) = 7 (0) = 5(0) = k (0) = s (0) = 0, f3 (0) = V2, h /i (0) = 1 preserves the initial data: 
ip (x, 0) = ^ (£, 0) and simplifies the solution flO)-fl3T7l). 



8. Two Soliton Solution 



Two well-known solutions of ( 17. ip are given by [78], |79j: 

AT 



and 



a cosh 3X + 3e 8iT cosh X 



^ 2 (X,T) = 4e u ! . (8.2) 

v ' ; cosh4X + 4cosh2X + 3cos8T v ; 

Use of the transformation (I7.12p . results in one and two soliton solutions for the nonautonomous 

nonlinear Schrodinger equation ( 12. 2ft . respectively. (See [78], [79] and jlQlj . for more details; N- 

soliton solutions are also discussed in [2], [62] and [7T].) 



9. Transformation of the Lax Pair and Zakharov-Shabat System 



If needed, an equivalent (nonisospectral) Lax pair for the nonautonomous Schrodinger equation 
( 12. 2p . which is discussed in Refs. [9], [H], [19], [21], [M], [85] for important special cases (see also [7] 
and [16]), can be derived, in general, from ( I6.3p -( }6T4"1) by inverting our transformation ( 17.121) (see [3], 
[55] for more details). In this paper, the required integrability condition ( 12 .3p (found for the soliton- 
like solution in Ref. [87] ) has been already incorporated into this transformation. Computational 
details are left to the reader. 



10. Examples 



A few simple examples are in order. (More examples can be found in [30], [87]; see also references 
therein.) 



NONLINEAR SCHRODINGER EQUATIONS 9 

10.1. Example 1. As noticed by Clarkson [21] (see also [2] and [68]), the equation 

^ + ^ = ^xV + 2H 2 ^ cu^O (10.1) 

does not pass the Painleve test and, therefore, is not integrable. The corresponding characteristic 
equation, ji" + u 2 fi = 0, has two standard solutions: 

2 

/i = — sin cut, = cos cut. (10.2) 

cu 

By our Lemma 1, a modified equation for harmonic solitons [45], [87] : 

#, + *. = Pi- + . m ft r ( °'f ' 0) : M V (10.3) 
4 4a (0) sin cut + cu cos cut 

(h , a (0) , /3 (0) 7^ and /i (0) are real constants) can be transformed into the standard form 
and hence is integrable. Here, \i (t) = \i (0) [4a (0) sin cut + cu cos cut] /cu and general solution of the 
corresponding Riccati-type system is given by 

cu4a (0) cos cut — cu sin cut ,„„ JS 

a(t) = H , 10.4 

w 4 4a (0) sin cut + cu cos cut v ' 

Pit) = „ " m , (10.5) 

4a (0) sin cut + cu cos cut 

2 (0) sin cut 

7 t = 7 , ; , 10.6 

' w IK J 4a (0) sin cut + cu cos cut' v ' 

5(t) = „ UH0) , (10.7) 

4a (0) sin cut + cu cos cut 

w w 4a (0) sin cot + co cos cot v y 

<5 2 (0) sin cut 

K (t) = K (0) , n . . ' . 10.9 

w w 4a (0) sin cot + cu cos cot v y 

Letting /i (0) = (5 (0) = 1 and a (0) =7 (0) = 5 (0) = e (0) = k (0) = 0, we arrive at the simple 
substitution: 

e -(^/4)x 2 tan^ / x tan CUt ^ 



V (z, t) = X -, , (10.10) 

V cos cut \ cos cut cu / 

which transforms (110. 3p into (12. 5D . (Replace cu = %k for the original equation in Ref. [43].) 
In a similar fashion, the following equation 

*, + + - «) * = (k+ ^XkL( 0) w * (iau) 

(A; 7^ 0, a (0) , /3 (0) 7^ and ^ (0) are real constants) is integrable. Indeed, the characteristic 
equation and the standard solutions are given by 

fi" + 4kfi' = 0; ^ = _L (1 - e ~«*) , Ml = l (10.12) 

and 

H = ^M- [A; + 2a(0) + (A;-2a(0))e- 4fc '] , A = e~ 2kt . (10.13) 
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The simplest case occurs when k + 2a (0) = [21J. Further details are left to the reader. 

10.2. Example 2. A soliton moving with acceleration in linearly inhomogeneous plasma was dis- 
covered in Refs. [18] and [19] (see also [9], [92]). For a modified equation with the integrability 
condition (12. 3p . 

,_ + _ + 2^ = 3-^ |^| ^ (10.14) 
where k, h , ao, (3 and /i are constants, we get /i (t) = /z (1 + 4a ^) and 

2/3 t (g + fct) /A 3 t (<5 + H) 2 

£ t =£ ° 1 i /i 7 ' K W = ^o-^ . — : , ■ (10.17) 

1 + 4a t 3 1 + 4a t 

The classical case [32], [H] corresponds to a = and h /j, (3l = —2 (with — > — fc). The simplest 
transformation: 

ij (x,t) = e < 2to - 4fc2 * 3 /3) x ( x _ 2k?, -t) , (10.18) 
when a = a = 7 = 5 = e = k = fi = and (3 = 1, is due to Tappert 



10.3. Example 3. The nonlinear Schrodinger equation [87J: 

■ dip d 2 ip g [3l h fi (3 2 2 hnim 

i— — — - — 7t = — o z W + : t^w\ (10.19) 

dt dx 2 (l + 4a t) 2 l + 4a t m Y ' V ' 

where 



/V + 2 ( 7o - (/3 2 - 4a 7o) V 



1 + 4a t 

with the help of the gauge transformation: 



(10.20) 



(1 + 4a t) 



^ = e -/W x ( x , t ) , = 2^2/ Z\ ,; 3 ' u/ , (10.21) 

can be transformed into a similar form 



% + X**-7TT-; — JT3* = i | 77^ Ixl X- (10.22) 



9oP 3 x = fop/Xpfl 
;i + 4a t) 3X ~ l + 4a t 
Looking for a traveling wave, we indicate the following soliton-like solution 

V |/u | (l + 4a i) V Hq \ J 



where 



s / xt \ = a ° x2 t Po xy t ( 7 o ~ W ~ 4c *o^o) f ) y 2 ( 1Q 24 ) 

1 + 4aot 

, «2, 2 7o - (£o - 8"o7o) * 

+90 ^ — — " 
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(a , P , 7 , S , e , k , fi , g , h , y are real constants) and the soliton profile is defined, as a 
solution of the second Painleve equation [89J, in terms of the nonlinear Airy function Ak Q (Q with 
asymptotics given by 

k Ai (C) , C -> +00 
A k>(0 = { .iAi-i/4 ci „( o( s\_n fir (10-25) 



r |C|- 1/4 sin ( a (C) - 9 ) + o (|C|" 1/4 ) , C -00. 



Here, Ai(£) is the Airy function, — 1 < k < 1 provided k 7^ 0, r 2 = — ix 1 In (1 — k^) , 

S (C) = ^|C| 3/2 -i 2 ln|C| (10.26) 

and 

9 = ^r 2 In 2 + arg T (l - l -r 2 ^j + ^ (1 - 2sign (A )) . (10.27) 

These asymptotics were found in [80] and had been proven rigorously in [3H], [39] (see [5], [TO] , 
[22], [23], [25], [27], [90] and references therein for study of this nonlinear Airy function). 

It is worth noting that, in a contrast to the previous case, this A-soliton moves with a constant 
velocity when a = (notice that the external field has essentially changed the soliton shape; see 
Ref. [87] for more details). 



11. Conclusion 



We have shown how to transform the nonautonomous and inhomogeneous nonlinear Schrodinger 
equation (12.2p ~ (j2.3p into the standard autonomous form (12. 5p that is completely integrable by the 
inverse scattering approach. This transformation is explicitly written in terms of Green's function 
of the corresponding linear problem (generalized harmonic oscillators); see Lemma 1 and equations 
(I3.ip - (l3.7p and f l3.8p - fl3.13p . Combination of these advances allows one to solve the Cauchy initial 
value problem for the generic nonautonomous integrable quantum system under consideration. 
Simple examples are considered and further examples can be found elsewhere (see, for instance, 
[30] . [87] and references therein). 
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